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O , ' The nonperturbative real-time evolution of quantum fields out of equilibrium is 

fl^ ' often solved using a mean-field or Hartree approximation or by applying effective 

r^ action methods. In order to investigate the validity of these truncations, we imple- 

ment similar methods in classical scalar field theory and compare the approximate 
dynamics with the full nonlinear evolution. Numerical results are shown for the 
early-time behaviour, the role of approximate fixed points, and thermalization. 
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>0 ■ 1 Testing truncations? 

^P \ An understanding of the nonperturbative evolution of quantum fields away 

from equilibrium is needed in many physical situations. Canonical examples 
are the universe at the end of inflation and the early stages of a heavy-ion 

p^ i collision. In both cases one would like to calculate how the energy initially 

contained in the inflaton or heavy ions is redistributed over the available de- 
grees of freedom, leading to a hot universe and a thermal quark-gluon plasma 
respectively. 

Jr~^' Because of the inherent real-time nature and the nonequilibrium setting, 

(-H I standard euclidean lattice methods are not applicable and approximations have 

to be introduced. In the last ten years or so a lot of attention has been 
given to mean-field approximations, in which the dynamics of a mean field 

S^ ■ coupled to Gaussian (quadratic) fluctuations is solved self-consistently. Well- 

JH \ known examples are the Hartree and leading-order large Nf approximations, 

with N f the number of scalar or fermion fields. Several paths going beyond 
homogeneous mean fields have been explored recently as well. For a list of 
references, see [1]. It goes without saying that the use of approximations 
will introduce deviations from the true evolution in quantum field theory. In 
order to gain trust in a particular method, it would be helpful to compare 
the approximate dynamics with the exact one. Unfortunately, in QFT such 
tests are not easy since the exact nonperturbative evolution is not known. 
Therefore we decided to use classical fields, regularized on a lattice, instead. 
In many aspects classical fields are similar to quantum fields, e.g. scattering is 
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present and the role of the thermodynamic hmit at fixed lattice spacing can be 
investigated. Furthermore, most approximation methods can be implemented 
as well, as they do not involve h directly. Last but not least, the 'exact' 
evolution can be calculated numerically which allows for a direct comparison. 

2 Classical Hartree and beyond 

We consider classical (p* theory in 1 + 1 dimensions. The dynamics is deter- 
mined by the equation of motion df(j){x,t) = {d^ — rn?)(j){x,t) — \(j)^{x,t)/2, 
given initial conditions c/'(x,0) and 7r(x,0) = dt4'{x,0) that are taken from a 
probability distribution p[7r, ^]. Ensemble averages are denoted with (•). The 
exact nonlinear evolution can be calculated by sampling initial conditions from 
p and solving the e.o.m. numerically for each of them. 

A classical Hartree- type of approximation can be implemented by replacing 
X(f)^ with 3X(j){<p'^) in the equation above. For translationally invariant systems 
{(f)'^) = {(j)'^{x,t)) is independent of x and the Hartree e.o.m. can be written 
conveniently in momentum space as d^(j){q,t) — —Lu'^(j){q,t). The eS'ective 
frequency squared is iUg = u}g + |A(0^), with w^ — q'^ +in'^. The unequal-time 
two-point function S{x — y; i, t') = (0(x, t)(j){y, t')) obeys in this approximation 
the usual mean-field equation of motion {df+u}i)S{q] t, t') = 0. For the Hartree 
approximation with inhomogeneous mean fields in the quantized version of this 
model, see [2]. 

Beyond the Hartree approximation unequal-time formulations typically 
become nonlocal in time. This is a disadvantage when numerical solutions are 
required. Therefore we reformulate the Hartree dynamics in terms of three 
basic equal-time two-point functions: G,p^>(x — y,t) = {xp{x,t)ip'iy,t)), with 
tjj = {tt, <j)}. In terms of these the Hartree evolution equations read 

dtG^^{q,t) = 2G^^{q,t), 

dtC^^iq, t) = -LU^G^^iq, t) + G^^{q, t), (1) 

dtG^^{q,t) = -2LJlG^^{q,t). 

These equations conserve a^^{q) = Gci,ci,{q,t)GT„T{q,t) — G^<{q,t) for each q. 

The evolution equations beyond Hartree include dynamical equal-time 
four-point functions. In the case of an 0{Nf) theory they contain all l/Nf 
corrections. Details can be found in [1]. It remains to specify the initial prob- 
ability distribution. We take a Gaussian ensemble with {ip{q,0)ip'{q' ,0)) = 
G^V' (9: 0)2TrS{q + q')5^^> , and 

G^^((Z,0)-ro/(g2+TO2), C(g,0)=To. (2) 

Note that this corresponds to the equilibrium ensemble when A = 0. 




Figure 1: Early-time evolution in the Hartree approximation (LO), with the inclusion of four- 
point functions (NLO), and 'exact' (MC). Left: mean field squared {(j)'"^). Right: zero-mode 
G'.Aq = 0,t). The initial temperature of the Gaussian ensemble is Tq = 5. 



3 Early time, fixed points, and thermalization 

Typical behaviour at early times is shown in Fig. |l| for two observables. We 
see that the LO evolution is in qualitative agreement with MC result. At NLO 
the agreement is impressive. For numerical details I refer again to [1]. Note 
that primes indicate dimensionless quantities and Tq = 3XTo/m^. 

An important issue in nonequilibrium field theory is the approach to equi- 
librium. One of the crucial tests for a truncated evolution is whether it is able 
to describe the thermalization regime or if the approximation breaks down be- 
fore. In a classical field theory, a convenient two-point function to investigate 
is G^7r(Q,i) = T{q,t). Using the equilibrium partition function as a guide, we 
see that this correlator can be interpreted as the effective temperature for a 
momentum mode q. In equilibrium, T{q) = T for all q. In Fig. (left) we show 
the time evolution of T(g, t) for three momentum modes (NLO, MC only). It is 
clear that the modes have a different effective temperature. The system seems 
to be quasi-stationary but is not thermal. 

A remarkable feature is that this behaviour can be understood completely 
in terms of fixed points in the truncated evolution equations. The stationary 
points in the Hartree approximation are easily found from Eq. (pi), and obey 



G:jq)^LO*/G;^iq) 
for the conserved combination at the fixed point 



g;^?) = 0, 



(3) 



with w*^ = Wg-|-|A(0^)*. These relations can be completed with the expression 



'-iq)=TQ'/u;'^^GUq)G;4q). 



(4) 
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Figure 2: Nonthermal behaviour. Left: effective temperature T'{q,t) for 3 momentum 
modes. Right: time-averaged temperature profile between < mt < 50, the fourth line is 
the analytic expression (a) at the fixed point in the Hartree approximation (Tq = 5). 



The first equality follows from the initial ensemble (0). Eqs. 
the fixed point completely. We find 
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(the second expression is a gap equation for 
of elliptic functions) and 



and can be written in terms 
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i.e., a nonthermal profile of G*^((7) at the fixed point. In order to see whether 
this fixed-point profile is realized in the dynamical evolution, we show the 
time average of G^Tr{q,t) for the Hartree, NLO, and full evolution in Fig. || 
(right). It is clear that the analytic expression not only describes the Hartree 
result, but also the profile from the NLO and in particular the full nonlinear 
evolution surprisingly well. This can be interpreted as a justification for the 
use of a Hartree approximation in this time interval. 

The fate of the fixed point can be determined by being patient. In Fig. ^ 
(left) we show the evolution of T{q = 0,t), for a higher initial temperature Tg 
as before. Also shown are the average temperatures T{t) = N^^'Y^ T{q^t). 
The Hartree evolution remains to be controlled by its fixed point. In the full 
nonlinear evolution on the other hand, all momentum modes obtain the same 
temperature. The change from a fixed-point profile to a thermal profile is 





Figure 3: Fate of the fixed point. Left: zero-mode temperature T'(q = 0,t) for LO and MC 
(Tg = 15). In the Hartree approximation the fixed point governs the dynamics for late times 
as well, whereas the full nonlinear dynamics shows thermalization. The average temperatures 
T'{t) appear as straight lines. Right: time-averaged snapshots of the temperature profile 
(MC only, T^ = 20, line 1: < mt < 1500, . . ., line 6: 13500 < mt < 15000). The fixed-point 
shape slowly disappears and the spectrum becomes thermal (flat). 

shown in Fig. ^ (right) for the fuh evolution only. Due to instabilities, it was 
not possible to reach such late times using the NLO evolution [1]. 

4 Summary 

In order to gain insight in how well truncated dynamics reproduces the full 
evolution in nonequilibrium quantum field theory, we studied the corresponding 
classical problem. We stress that the aim was not to use classical fields to 
approximate the quantum theory, but rather to study truncated dynamics in a 
field theory where the exact evolution is available. Our findings are that while 
at early times the truncated evolution works remarkably well and the outcome 
can be characterized nicely in terms of fixed points of the associated Hartree 
equations, it breaks down before the thermalization regime is reached. 
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